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ABSTRACT 

The global star formation rate has decreased significantly since z ~ 1, for reasons 
that are not well understood. Red-sequence galaxies, dominating in galaxy clusters, 
represent the population that have had their star formation shut off, and may therefore 
be the key to this problem. In this work, we select 127 rich galaxy clusters at 0.17 < z < 
0.36, from 119 square degrees of the Canada- France-Hawaii Telescope Legacy Survey 
(CFHTLS) optical imaging data, and construct the r'— band red-sequence luminosity 
functions (LFs). We show that the faint end of the LF is very sensitive to how red- 
sequence galaxies are selected, and an optimal way to minimise the contamination from 
the blue cloud is to mirror galaxies on the redder side of the colour-magnitude relation 
(CMR). The LFs of our sample have a significant inflexion centred at M r > ~ —18.5, 
suggesting a mixture of two populations. Combining our survey with low redshift 
samples constructed from the Sloan Digital Sky Survey, we show that there is no 
strong evolution of the faint end of the LF (or the red-sequence dwarf-to-giant ratio) 
over the redshift range 0.2 < z < 0.4, but from z ~ 0.2 to z ~ the relative number 
of red-sequence dwarf galaxies has increased by a factor of ~ 3, implying a significant 
build-up of the faint end of the cluster red-sequence over the last 2.5 Gyr. 
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It has been established that the galaxy population today 
can be broadly divided into two categories: those that have 
red colours, consisting of mostly non-star- forming galax- 
ies (the 'red-sequence'); and those with blue colours and 
active star formation (the 'blue cloud'). This colour bi- 
modality c an be modelled by a s um of two Gaussia n dis- 
tributions <|Balogh et al.l l2004al lbt baldrv et all I2004l 1. and 
it has been observed out to z ~ 1 ( Bell et al J 12004V The 
red population exhibits a tight correlation between colour 
and magnitude, with brighter g alaxies being r edder. Due to 
the age-metallicity degeneracy (|Worthevlll994l ). the slope of 
this colour-magnitude relation (CMR) of the red population 
could be attrib uted to either an age or met allicity sequence. 
The study by iKodama fc Arimotol l| 19971 ) concluded that 
metallicity variation dominates the slope, because of its rela- 
tively slow evolution. Nonetheless, some age variation along 
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the red-sequence is also detected (e.g. Trager et al. _ 2000l ; 
iNelan et alj|2005l : ISmith et al.ll2006l : lAllanson et al.l l2009). 



The origin o f red-sequence gal axies i s still an ope n ques - 
tion. Studies by iBell et all (|2004) and iFaber et all (|2007t ) 
concluded that brighter red-sequence galaxies are built up 
through dry mergers since z ~ 1, based on observations that 
the B— band luminosity density of red-sequence galaxies re- 
mains constant since z ~ 0.9. With the dimming of galaxies, 
the luminosity density would be overproduced if the num- 
ber density of red-sequence galaxies remained constant, as 
in the pure passive evolution scenario. H owe ver, other stud- 
ies by , for example. ICimatti et alj (2006) and lScarlata et all 
(2007), showed that the number density of the massive red 
galaxies remains constant and their luminosity function is 
consistent with passive evolution. It is the number density of 
less-massive galaxies that decreases rapidly with increasing 
redshift, and this can be explained by a gradual quenching 
of star formation in those galaxies. Therefore, no significant 
amount of dry mergers are required to explain the formation 
of massive red-sequence galaxies. Which of the two proposed 
scenarios is more important is still not clear. 
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As to how the faint end of the red-sequence builds up, 
it is even less clear. Especially in galaxy clusters, different 
studies have yielded conflict i ng re s ults. For example , ob- 
servations bvlDe Lucia et alj (l2007T ). [stott et all (|2007t) and 
iGilbank et al'l |2008f) found a significant deficit of faint red- 
sequence galaxies in high redshift clusters, compared to low 
redshift clusters. This supports a scenario where low mass 
galaxies have their star formation shut off and move onto 
the red-sequence recently (z f). However. lAndreonl (120081 ) 
and lCrawford et alj (120091 ) found no evolution of the faint- 
end slope spanning < z < 1 .3, which implies the build-up 
of the faint red-sequence galaxies was complete d by z ~ 1.3. 

Current galaxy formation models (e.g. iBower et alj 

2006) produce too many red galaxies in c l usters and groups 
that a re not observed JWolf et alj 120051 ; IWeinmann et all 
l2006al lbl; lBaldrv et al.ll2006h . indicating a fundamental prob- 
lem in our understanding of how star formation is quenched 
in dense environments. More precise observations are needed 
to resolve this problem. 

In this work, we make use of optical imaging data 
from the Canada-France-Hawaii Telescope Legacy Survey 
(CFHTLS) data to construct a large sample of 127 clusters 
and study their red-sequence luminosity functions over the 
redshift range 0.17 < z < 0.36. In SjHand 33 we describe the 
data, and the cluster detection algorithm. The properties 
of our cluster catalogue, and a local comparison sample are 
described in 3Hand 321 We present our methods for measur- 
ing the red-sequence luminosity function and dwarf-to-giant 
ratio in 3H1 and the results in 33 ln 3H we discuss the compar- 
ison with literature and possible systematics, and conclude 
in 3U 

Throughout this paper, we assume a cosmology with 
O m = 0.3, Oa = 0.7 and H o =70 km s" 1 Mpc" 1 . All magni- 
tudes are in the AB system unless otherwise specified. 



Table 1. Total area covered by the Wide survey and area used 
in this work. 



2 DATA 

2.1 The Survey 

The CFHTLS is a joint Canadian and French imaging sur- 
vey in it*, g' , r', i' and z filters using MegaCam, with an 
approximately lxl square degree field of view. We use the 
"Wide" and "Deep" surveys in this work. The Wide survey 
covers a total area of 171 square degrees with a total expo- 
sure time of about 2500 seconds in g' band, 1000 seconds in 
r' band and 4300 seconds in i' band per pointing. The Deep 
survey contains four 1 square degree fields with exposure 
times ranging from 33 hours in u* band to 132 hours in i' 
band. 

The survey started in 2003 and is now complete, but 
not all data have been processed and released yet. The data 
we used in this work are from data release T0004 (released 
internally July 3, 2007). We carry out our cluster detection 
using the Wide survey data over a total area of 119 square 
degrees (details summarized in Table [TJ, where the photom- 
etry in all three of the filters, g' , r' and i' , is available. We 
use data from one of the four Deep fields, Dl, as auxiliary 
data to examine surface brightness selection effects ( j)2.2.1|) 
and other possible systematics ( i]8.3.2[) . 



Wl 


W2 


W3 


W4 


Total (sq. deg) 72 


25 


49 


25 


Used in this work (sq. deg) 42 


20 


41 


16 



2.2 Photometry 

For our analysis, we are interested in measuring the colours 
and total magnitudes of galaxies. The photometric cata- 
logue we use was produced by TERAPIX (Traitement El- 
ementaire, Reduction et Analyse des PIXels de megacam). 
Magnitudes in the catalo gue are measured using SExtractor 
ijBertin fc Arnou"tslll996t ) within different apertures. We take 
magjoLuto, where flux is measured within the Kron radius, 
as the total magnitude of galaxies. 

We measure colours of galaxies using magnitudes within 
a fixed aperture of 4.7 arcsec in diameter. The seeing*! of 
all images in this release is better than 1.3 arcsec, and the 
typical seeing of each pointing is ~0.94 arcsec in g' , ~0.86 
arcsec in r' and ~0.81 arcsec in i'. For each pointing, the 
maximum difference in the seeing of the stacked, single-filter 
images between g' , r' and i' band is less than 0.5 arcsec in 
the worst case — small compared with the radius of the 
aperture used to measure the colour. Therefore, we do not 
convolve the images in different filters to the same seeing. 
TERAPIX compared the photometry of objects in the over- 
lap regions bet ween the CFHTLS and the Sloan Digital Sky 
Survey (SDSS) (|York et alJl200Ch . and found that the mean 
offset in g , r',and i' filters in each individual pointing ranges 
from ~ 0.01 to ~ 0.05 mag. 

The galac t ic E( B-V) extinction calculated using the 
ISchleeel et all fl998) dust map is provided in the cata- 
lo gue for each object. We use th e A/E(B-V) valu es given 
in lSchlegel et~ai] fl99Sl ) in SDSS |York et al.ll200ol ) it, g, r, 
i, and z filters to calculate the extinction correction and ap- 
ply this to the observed magnitude. The median extinction 
in r' band is 0.072 mag in Wl, 0.056 in W2, 0.032 in W3 
and 0.217 in W4. 



2.2.1 Depth and Photometric Uncertainty 

To examine how surface brightness might affect the com- 
pleteness of the sample, we plot the surface brightness for 
the brightest pixel, fimax (which is a good proxy for the 
central surface brightness, measured using SExtractor), vs. 
magnitude for galaxies that are in the overlap region be- 
tween Dl and Wl, as shown in the left column in Figure [1] 
Blue points represent objects in Dl and red points Wl. In 
g' band, at about g' ~ 25, Wl starts to become incomplete 
due to the surface brightness detection limit. To quantify 
this, we calculate the fraction of sources in Dl that are also 
detected in Wl as a function of magnitude. As shown in the 
right column of Figure [T] Wl is ~ 100 per cent complete 
at g — 25, and r' = 24. Since the (</ — r') colour of red 



* denned by TERAPIX as twice the median half-light radius of 
a selection of point sources on each CCD as measured by SEx- 
tractor. 
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Figure 1. Left column: central surface brightness, 

P^vaax , as a 

function of g' (top) and r 1 (bottom) magnitude. Red points are 
objects in Wl and blue points Dl. Right column: fraction of 
sources in Dl that are also detected in Wl, as a function of mag- 
nitude. 



galaxies at the highest redshift we focus on in this paper is 
around 1.8 (Figure 0), we make a cut at r' = 23.2 to ensure 
completeness in the g' band. 

To determine the uncertainty on the colour, we examine 
the difference in the colours measured from each pointing, 
for galaxies in the overlap regions between pointings. "fl In 
Figure [2] we plot the standard deviation of (g' — r') (left 
panel) and (r' — i') (right panel) colour differences as a 
function of magnitude, and fit the logarithm of it by two 
straight lines. At the magnitude limit of r' = 23.2 of our 
catalogue, the uncertainty on the (g 1 — r') colour is about 
0.2 mag, which makes isolating red-sequence galaxies from 
the contamination of the blue cloud less reliable at fainter 
magnitudes (see [JS]). 



2.3 Star-Galaxy Separation 

TERAPIX classified objects in the magnitude range 17.5 < 
i 1 < 21.0, using the stellar locus in the half-light radius vs. 
magnitude plot, in a series of 10 arcmin cells distributed 
over each MegaCam stack. All objects fainter than i' — 
21.0 are considered galaxies. At magnitudes brighter than 
i! — 17.5, we keep objects that have SExtractor parameter 
class_star < 0.98, but exclude those on the stellar locus, as 
galaxies. Objects in our final galaxy catalogue are plotted as 



0.01 



Figure 2. Left panel: the standard deviation of (g' — r') colour 
differences in overlapping pointings as a function of r' magnitude, 
fitted by two straight lines on logarithmic scale. Right panel: the 
standard deviation of (r' — i') colour differences as a function of 
i' magnitude. 
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Figure 3. Half-light radius vs. magnitude. Points represent a 
subset of the galaxies that are in our final galaxy catalogue, and 
circles represent stars. See text for details on star-galaxy separa- 
tion. 



small points in Figure|3] and stars are represented by circles. 
Only a subset of the data is plotted for clarity. 



tThe residual mean offset in colour between two overlapping 
pointings is < 0.05 mag, and we correct for this residual offset 
before we stack all the overlap regions together to examine the 
distribution of the difference between colours measured from two 
overlapping pointings. 



3 CLUSTER DETECTION 

The cluster detection method we use is ba sed on the spirit 
of the cluster-red-sequence (CRS) method l|Gladders fc Yed 
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2000), but with some modifications. The foundation of the 
CRS technique is that in every cluster there is a popula- 
tion of red galaxies that forms a tight red-sequence in the 
colour- magnitude diagram (CMD), which can be used as an 
overdensity indicator. Four major steps are involved in the 
cluster detection: 1) model the CMD at different redshifts; 
2) select a subsample of galaxies that belongs to each colour 
slice; 3) count galaxies around a grid of positions in the sky 
and estimate the background in the same way; 4) select over- 
densities as cluster candidates. Below, we describe each step 
in detail. 

3.1 Model CMD 

First we model the CMD for the red-sequence galaxies. Since 
we only need a model that produces the correct CMR over 
a small redshift range, and we will calibrate it empirically, 
a si mple model is suffici e nt. W e use a single burst model 
from iBruzual fc Chariot! l|2003h with a Salpeter IMF and 
formation redshift of 2/ ~2. We follow the passive evolu- 
tion of galaxies with different metallicities, and calibrate the 
model by defining a magnitude-metallicity relation so that 
it reproduces the observed CMD o f E+SO galaxies in Coma 
(slope -0.0743, zero point 4.21, iBower et ail 1 19921 ). The 
resulting magnitude-metallicity relation is in good agree- 
ment with the result from a more complicated model by 
iKodama fc Arimotol (|l997h . 

Under the assumption that this magnitude-metallicity 
relation does not evolve with redshift, we can define for each 
metallicity a correlation between the model colour and mag- 
nitude for the passively-evolving galaxies at each redshift. 
We fit the resulting CMR with a straight line. To model m* 
at different redshifts, w e take the m* at z — estimated 
bv lBlanton et al.l l|200ll) using SDSS data, transform it into 
CFHTLS filters, and evolve it using the same model with 
solar metallicity. 

Figure|4]shows the CFHTLS filter transmission curves^, 
and into which filters the 4000A break (the most prominent 
feature in a red-sequence galaxy's spectrum) falls at differ- 
ent redshifts. When the 4000A break is in the g' filter, the 
(9' ~ r ') colour is very sensitive to a small shift of the break 
position, and thus gives the best redshift resolution. When 
the break is approaching the boundary between g' and r' 
filters, the (g' — r') colour becomes degenerate. Therefore, 
we limit our study in this work to z < 0.36. 

Figure [5] shows the modelled (<?' — r') vs. r' diagram 
as a function of redshift, with m* indicated. The redshift 
interval between slices is 0.03. The m* and colour at m* 
are tabulated in Table [2] with a small, empirical correction 
applied to the redshifts as discussed in $4.3\ 



3.2 Subsamples in each Colour Slice 



The next step is to construct a subsample of galaxies for each 
redshift based on the colours of the galaxies. At each red- 
shift, we define a colour slice centred on the model CMR, 
with the half width determined by the in trinsic scatter in 
the c olours of red-sequence galaxies, 0.075 l|Gladders fc Yed 
2000). The colour slices overlap with each other to avoid 
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Figure 4. Position of the 4000A break at different redshifts in 
CFHTLS filters. Solid curves are total filter transmissions (fil- 
ter+mirror+optics+CCD), and triangles indicate the positions 
of the 4000A break at different redshifts. 




Figure 5. Model CMD in CFHTLS filters as a function of redshift 
with m* indicated. The redshift i nterval is Az = 0. 03. The z ~ 
relation is calibrated to Coma llBower et al.lll992h . The CMD 
is evolved passively, with a fixed magnitude-metallicity relation. 
Points represent a subset of galaxies from Wl that are associated 
with the colour slice at z = 0.25 (see £|3.2|l . 



missing clusters on the edge of each colour slice. We then 
calculate, for each galaxy, the probability that it belongs to 
a specific colour slice, assuming that the probability distri- 
bution of the colour of each galaxy is Gaussian, with a = Ac, 
dominated by the colour errors estimated from the overlap- 
ping regions as shown in Figure [2] For each colour slice, all 
lhttp://www3.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/megapipe/docs/filters.^feiliabdes with probabilities larger than 10 per cent are se- 
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Table 2. Model m* and colour at m* as a function of redshift. 



z 


m*, 


(g' —r') at m*, 


0.19 


17.65 


1.048 


0.22 


18.1 


1.128 


0.25 


18.48 


1.213 


0.28 


18.81 


1.309 


0.31 


19.12 


1.394 


0.34 


19.4 


1.462 


0.37 


19.67 


1.517 



lected. A subset of the galaxies from Wl that are associated 
with the colour slice at z = 0.25 are shown as points in 
Figure [5] At the faint end, due to the larger colour uncer- 
tainties, galaxies belonging to that colour slice spread out of 
the boundary of the slice. 

This colour-weight ing step is the same as prescribed by 
iGladders fc Yeel (|2000r ). In the original CRS technique, af- 
ter the colour weighting, a magnitude weighting is applied 
to downweight the fainter galaxies because the contrast be- 
tween cluster and the field is lower at faint magnitudes. We 
do not want to bias our detected clusters to a certain lumi- 
nosity function shape, so we do not apply any magnitude- 
weighting here. However, to avoid high contamination from 
the field at the very faint magnitudes, we only use galaxies 
brighter than m* + 2 for detection. For the purpose of our 
study here, we only need a sample of rich clusters without 
worrying about the completeness of poorer systems, which 
further justifies the modification we make here. 



3.3 Significance Map and Detection 

In each colour slice, we count the number of galaxies in the 
subsample within a circle of radius ~ 0.5 Mpc (the typi- 
cal size of a cluster core) around a grid of positions ~ 0.1 
Mpc apart in the sky. This gives the cluster+field count 
at each position. The field count is estimated in the same 
way, but from the average of 500 random positions in the 
sky. The difference between the total number count and the 
background count at each grid position divided by the rms 
deviation of the background distribution of the 500 random 
positions (comparable to Poisson uncertainty), 07, indicates 
the significance of the overdensity smoothed on a 0.5 Mpc 
scale at that position, i.e. 

U= ^ '- . (1) 

FITS images of the significance maps for each colour 
slice are created. SExtractor is then run on those images to 
pick out peaks on the significance maps. We keep everything 
that has at least 1 pixel above 5cr as cluster candidates. 

We now have a list of crude positions and significance of 
cluster candidates. If multiple peaks are detected within ~1 
Mpc from each other, either in the same colour slice or ad- 
jacent slices, only the one that has the highest significance 
will be kept. In the following sections, we refine the clus- 
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Figure 6. Evolutionary tracks of (</ — r') vs (r 1 — i') with red- 
shift for three different p opulations produced by the model of 
iBruzual fe C hariot (2003). Squares are a population that formed 
10.3 Gyr ago and have constant star formation. Pentagons rep- 
resents a 8 Gyr old single-burst model and triangles a 10.3 Gyr 
single-burst model. The lines are to guide the eyes. 

ter catalogue, by using additional colour information, and 
improving the centring and redshift resolution. 

3.4 Refinements to the Cluster Catalogue 

3.4-1 High-z Contamination 

For any cluster candidate we detected, higher-redshift blue 
galaxies projected along the line-of-sight may have the same 
(9' ~~ r ') colour as the red-sequence members of that clus- 
ter. This is demonstrated in Figure [5] which shows the 
evolutionary tracks of (g' — r') vs. (r' — i') with redshift 
for three different popula tions produced by the model of 
IBruzual fc Charlotl (|2003T ). Squares are a population that 
formed 10.3 Gyr ago and have constant star formation. Pen- 
tagons represents a 8 Gyr old single-burst model and trian- 
gles a 10.3 Gyr single-burst model. The lines are to guide 
the eyes. For the two old models, galaxies in the redshift 
range 0.4 < z < 0.9 all have similar (g' — r') colour. How- 
ever, they have different (r' — i') colours. Therefore, (r' — i') 
colour can effectively eliminate high redshift galaxies from 
the red-sequence galaxies at a lower redshift. This is further 
demonstrated in Figure [7] which shows the stacked CMDs 
of galaxies within 0.5 Mpc around each cluster in (</ — r') 
vs. r' (black points in the left column) and (r' — i') vs. i' 
(black points in the right column) planes at three redshifts. 
Pentagons indicate the position of m* . Central solid lines in 
the (g' — r') vs. r' plots indicate the model CMR and the up- 
per and lower bounds indicate the width of the colour slice, 
±0.075 mag. On the (r' — i') vs. i' plot, the central lines are 
the resulting CMR from the same model that produces the 
(9' ~ r ') vs - r ' relation. The half width of the colour slice 
in [r' — i') vs. i' is 0.2 mag (justified below) as indicated 
by the upper and lower solid lines. Blue points are galaxies 
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Figure 7. Stacked CMDs of galaxies within 0.5 Mpc around each cluster in (g' — r') vs. r' (black points in the left column) and (r 1 —i') 
vs. i' (black points in the right column) planes at three different redshifts. Pentagons indicate the position of m* . Central solid lines in 
the (g' — r') vs. r 1 plots indicate the model CMR and the upper and lower bounds indicate the width of the colour slice, ±0.075 mag. 
On the [r 1 — i') vs. i' plot, the central lines are the resulting CMR from the same model that produces the (g' — r') vs. r' relation. The 
half width of the colour slice in (r' — i') vs. i' is 0.2 mag as indicated by the upper and lower solid lines (see text for explanation). Blue 
points arc galaxies that belong to the red-sequence in (</ — r') vs. r' . As we can tell, many of them do not fall on the red-sequence in 
the (r' — j') vs. i' plane, which shows the advantage of using two-colour information to select red galaxies. 



that belong to the colour slice in the (g' — r') vs. r' plane. In 
the (r 1 — i') vs. i' plane, most of those galaxies still fall on 
the red-sequence at low redshift, but as we go to higher red- 
shift, many of them fall off the red-sequence even at bright 
magnitudes. Therefore, only galaxies that are on the red- 
sequence in both ((?' — r') vs. r' and (r' — i') vs. i' planes 



are considered as potential red-sequence cluster members. 
Note we only use (r' — £') colour to eliminate obvious back- 
ground contamination; therefore we take a wide colour slice 
in (r' — i') vs. i' so that we do not lose red-sequence galaxies 
due to the uncertainty on the (r' —i') colour and any slight 
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mismatches between the model (<?' — r' ) and (r' — i') c olour 
slices. This technique has been used in Andreor] (2006). 



3.4.2 Centres 

As described above, the centre of each cluster is determined 
by the position at which the number of red-sequence galax- 
ies within a radius of 0.5 Mpc selected using (g' — r') vs. r' is 
maximised. We refine the determination of the centre in sev- 
eral ways. We divide the region within a radius of 0.5 Mpc 
around each cluster further into smaller grids of size ~ 0.06 
Mpc and count galaxies that belong to the red-sequence in 
both (<?' — r') vs. r' and (r' — i') vs. i' planes in a circle of ra- 
dius ~ 0.1 Mpc. We also calculate the luminosity-weighted 
centre using galaxies that belong to the red-sequence in both 
(g' — r') vs. r' and (r' — i') vs. i' planes in a circle of radius 
~ 0.5 Mpc around each cluster. We compare the centres ob- 
tained in both ways with the position of the brightest cluster 
galaxy in each cluster. For some rich, symmetric systems, the 
three centres agree with each other well. In the calculations 
later in the paper, we use the luminosity-weighted centre. 



3.4.3 Redshift 

As described in i]3.1l fFigure[5"Tl. the redshift interval between 
adjacent colour slices is 0.03. To refine the redshift estimate 
of each cluster candidate, we insert two more colour slices 
between two existing adjacent ones. To determine, for each 
cluster, which redshift the model CMD fits the observed 
red-sequence the best, two criteria are used: 1) we count the 
number of red-sequence galaxies within a circle of 0.5 Mpc in 
radius around the luminosity-weighted centre of that clus- 
ter; 2) we calculate the deviation of the colours from the 
model CMR for all red-sequence galaxies that are brighter 
than m* + 2. We do this for each cluster in several adjacent 
colour slices and determine which one gives the highest num- 
ber count and least deviation in colour. If these two criteria 
give the same optimal colour slice, then the redshift of that 
colour slice is assigned to the cluster. If these two select two 
different optimal colour slices, the one that gives the least 
deviation is chosen only if at least five galaxies are used in 
the fit and the deviation is significantly smaller than that 
from the slice that gives the highest count. We check how 
well this works by stacking all clusters at the same redshift 
and plotting the observed (g' — r') vs. r' relation against the 
model. 

Finally, we interpolate the CMD in the (r' — i') vs. i' 
plane for the corresponding new redshift. A new centre is 
calculated for each cluster using the galaxies that fall on the 
red-sequence defined by the finely-interpolated colour slices 
in both the (g' — r') vs. r' and (r' — i') vs. i' plane. The 
significance of each cluster is re-estimated around the new 
centre. 



4 CLUSTER PROPERTIES 
4.1 Richness 

Since we used the number of red-sequence galaxies brighter 
than m*+2, and within only 0.5 Mpc from the cluster centre, 
in the cluster detection to reduce the noise, throughout this 
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Figure 8. Correlation between our N re d m *_|_2 measured within 
0.5 Mpc and the n umber of cluster me mbers within r2oo from the 
MaxBCG sample llKoester et alj|2007t) . 



paper we use this number, denoted as N re d, m *+2, to indi- 
cate the richness of the clusters in our sample. However, the 
typical extent of a cluster, T2oo (the radius within which the 
mean density is 200 times the critical density), is larger than 
0.5 Mpc; therefore, to get some idea of how our N re d,m'+2 
corresponds to the more commonly used richness indicator, 
A200 (the number of cluster members within r2oo), we plot 
N re d,m*+2 vs. N200 for the ten clusters that are both in our 
sample and the MaxBCG catalogue l|Koester et al.l 120071 ) 
at this redshift range in Figure [3 Based on the relation 

^200 



0.26A 2 °6o 2 Mpc (|Johnston et al.l 120071 ; lHansen et ail 
2007), the r2oo for clusters that are in the lower left corner 
of Figure [8] is about 0.8 Mpc; thus, their N re d,m*+2 as mea- 
sured within 0.5 Mpc is only slightly lower than A^oo- As one 
goes to richer clusters, the difference between N re d,m*+2 and 
A200 becomes larger. Nonetheless, we can approximately 
scale N red , m .'+2 to A200 using Figure 

The richest cluster in our catalogue has a N re d,m*+2 of 
46, a known Abel l cluster (Abell0362) with Richness Class 
1 and z = 0.1843 (|Cruddace et alj|2002l ). Its CMD is shown 
in Figure [9] Solid lines indicate the colour slice it belongs to 
with the star indicating the position of m* , and crosses are 
galaxies within a radius of 0.5 Mpc from the centre (before 
background subtraction). 



4.2 Mass Estimates 

To get some idea of how massive the clusters in our sample 
are, we compare the surface density of the clusters in our cat- 
alogue with that from the Hubble Volu me light cone. Base d 
on the MS (spheres) cluster catalogue (|Evrard et al.| [2002). 
there are 1.26 clusters per square degree with mass greater 
than 1.4 x 10 14 Af Q in the redshift range of 0.17 < 2 < 0.36, 
and this surface density corresponds to that of the clusters 
with N re d,m*+2 > 12 in our catalogue (as can be seen in 
Figure [TU)) . 
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Abell 0362 (z 3pee =0. 1843) 
= 0.20 




Figu re 9. The CMP of the richest cluster in our catalogue (Abell 
0362 ICruddace e t al. 2002) before background subtraction. Solid 
lines indicate the colour slice it belongs to with the star indicating 
the position of m* , and crosses are galaxies within a radius of 0.5 
Mpc from the centre. 
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Figure 10. Cumulative surface density of the clusters we de- 
tected as a function of N red m *+2 m the redshift range 0.17 < 
z < 0.36, in all four wide fields. Each line style represents one of 
the four wide fields. 



We have calculated the projected correlation function, 
uj(9), of the clusters in each of the four wide f ields sepa- 
rately, following the method of lLandv fc Szalavl (|l99of ). We 
select all clusters in our sample with 0.17 < z < 0.36, and 
a richness N re d,m*+2 > 10, which leaves us with 74, 31, 57 
and 32 clusters in the fields Wl through W4, respectively. 
This richness limit is slightly poorer than that adopted for 
most of the analysis in this paper in order to have sufficient 



statistics for the following analysis. This means our results 
in the paper actually correspond to systems slightly more 
massive than the limits given here. 

The correlation function is shown in Figure [TT1 The ran- 
dom point distribution we use for comparison does not ac- 
count for the intrinsic size of our detection filter; this means 
we will not measure the correlation function accurately on 
angular sizes of about 9 < 0.1 degree. We therefore fit the 
data with a power law function co(0) — A^9 1 ~~ 1 over the 
range 0.2° < 6 < 2° . We do not have enough data to mea- 
sure both the amplitude an d slope, so we fix 7 = 2.15 as 
measured locally for clusters (jGonzalez et alj|2002l ) . We find 
the best fit amp litude by maximizi ng the likelihood func- 
tion as given in lGonzalez et al.l (|2002r i. and take the 95 per 
cent confidence limits to be the value where the relative like- 
lihood is 0.1. These values are given in each panel in Figure 
1111 There is considerable variation from field to field, to be 
expected since the areas are still fairly small. The two fields 
with the most data, Wl and W3, actually represent the 
least- and most-clustered data, respectively, though we note 
that only the W3 data actually puts useful constraints on 
the clustering amplitude. 

Deprojecting the angular correlation function using the 
Limber equation and assuming our standard cosmology, and 
including a Gaussian fit to the observed N(z) distribution 
(see Figure I13[) . we find that the amplitude corresponds 
to a physical correlation length of r ~ 41 Mpc for W3, 
and 23.8 Mpc in the least-clustered field, Wl. This range 
of correlation lengths in turn corresponds t o a space den- 
sity o f n c = (0.2 - 4.6) x 10" 6 (Mpc) 3 (e.g. IColberg et al 



2000). We use the Millennium simulation (|Springel et al 
20051 ) output at z ~ 0.3, and find that this space density 
corresponds to clusters with mass M > 1.0 x 1O 14 M (Wl) 
and M > 4.1 x 10 14 Af e (W3). 

As an alternative to deprojecting the angular cor- 
relation function, we also compare our correlation func- 
tion directly with that of the Hubble volume simulation 
l|Evrard et al.l[20 02), by selecting clusters in the full-sky light 
cone, within a similar redshift range as the data and in dis- 
crete bins of total mass. We then compute the correlation 
function in the same way as for the data; the results are 
shown in Figure ITTI as triangles (M > 8.0 x 10 13 Mq), crosses 
(Af > 3.0 x 1O 14 M ) and squares (M > 5.0 x 1O 14 M ). 
In all fields except Wl, the comparison with the data sug- 
gests our sample with N re d,m*+2 > 10 is limited at about 
M > 2 x 10 14 Mq, consistent with the deprojection analysis 
above. 

The CFHTLS Wl field also overlaps with the XMM- 
LSS field and we recovered, in the redshift range we are 
interested in, all 9 X-ray confirm ed clusters in the XMM- 
LSS archive (|Pacaud et all l2007h . These clusters are all 
low-richness clusters in our sample. Their X-ray tempera- 
tures range from 1 keV to 3 keV. From the o bserved mass- 
temperature relation of g alaxy clusters (e.g. lEvrard et all 
119961 : IShimizu et all 120031 ) . the mass of a cluster with a X- 
ray temperature of ~ 2 keV is ~ 1 x 1O 14 M . This again 
indicates that the clusters in our sample have masses of 
M > 1 x 1O 14 M . 

Throughout this paper, unless otherwise specified, the 
analysis is carried out using the subset of 127 clusters 
that satisfy N re d,m*+2 > 12. We emphasize that such a 
richness-selected sample does not correspond exactly to a 
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Figure 11. For each of the four CFHTLS survey fields, we show 
the angular correlation function for clusters with 2Vj.edm*+2 > 
10, as the solid points with Poisson error bars. The solid line is 
the best fit power law w(0) = AuB 1 -^, with 7 = 2.15 fixed. The 
dotted lines show the 95 per cent confidence limits. The trian- 
gles, crosses and open squares are theoretical angular correlation 
functions for clusters with M = 8.0 X 1O 13 M , M = 3 X 1O 14 M , 
and M = 5 X 10 14 Mq, respectively. These were computed from 
the Hubble volume simulation over a similar redshift range as the 
data. 

mass-selected one, and thus a mass limit cannot be pre- 
cisely denned. Nonetheless the analyses above show that 
our results are applicable to fairly massive clusters, with 
M > 1 X 10 14 Mq; the sample is not dominated by low-mass 
groups. 



4.3 Redshift Accuracy 

The redshifts we initially assign to our clusters are based on 
the photometric data. In this section, we assess how good 
our photometric redshift is. 

As mentioned in £ 14.11 there are ten co mmon clusters 
in ou r catalogue and the MaxBCG catalogue |Koester et al.l 
2007). For those clusters, we compare the spectroscopic red- 
shifts of those BCGs with the photometric redshifts we as- 
signed to those clusters. This comparison indicates that our 
photometric redshift is systematically lower, by a constant 
amount of 0.0435. Therefore, we apply this constant shift to 
the estimated photometric redshifts from our initial model. 
(Recalib ration in this way is a standard part of the CRS 
method, iGladders fc Yeell200Ch All redshifts quoted in this 
paper have been adjusted in this way. Squares in Figure [121 
show the comparison between the adjusted z v hot and the 
Zs P ec of MaxBCG clusters. 

We also compare the adjusted photometric redshifts of 
the X-ray confirmed clusters fi ]4.2[) in our catalogue with 
their spectroscopic redshifts, shown as crosses in Figure [T2l 
The two agree with each other well, providing an indepen- 
dent verification of our correction. The rms of the difference 
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Figure 12. Comparison of our estimated photo metric redshift 
with the spectroscopic red shift from the MaxB CG jKoester et al.l 
120071) and X-ray sample llPacaud et al,l [20071) . Squares are our 
z phot vs - Zspec from maxBCG sample, with Zphot corrected by a 
constant shift of 0.0435. Crosses are our corrected Zp^ot vs. the 
Zspec of the X-ray confirmed clusters from XMM-LSS in Wl field. 
The two agree with each other well, providing an independent 
verification of our correction. 

between the z v y lotz and z spec (including both the MaxBCG 
and X-ray samples) is ~0.014, with very little bias. 

Figure [T3] shows the distribution of clusters as a func- 
tion of the adjusted z v hot in our sample, split into two 
subsamples based on N re d,m*+2- The solid histogram rep- 
resents the distribution of the subsample of clusters with 
N re d,m*+2 > 20, and the dotted histogram represents those 
with 12 < N re d,m*+2 < 20. The distribution is smooth with 
redshift, peaking at around z ~ 0.3. 



5 LOCAL COMPARISON SAMPLE 

It is useful to extend our redshift baseline by examining 
a sample of local clusters, using our methods consistently. 
Therefore, we make use of the low-redshift galaxy group cat- 
alogue by Yang et al. (2007, hereafter Yang07), selected us- 
ing a Friend s-of-Friends algor ithm from the SDSS spectro- 
scopic data |York et al.ll2000l ). We calculate N re d,m*+2 for 
these groups/clusters the same way as for our own cluster 
sample, except that the red-sequence galaxies are selected 
in (it — r) colour to bracket the 4000A break, because of 
the lower redshift of the sample. We select a subset of clus- 
ters that have N re d,m*+2 > 12, the same richness cut of our 
cluster sample. In order to carry out the calculations (e.g. 
background subtraction) exactly the same way as for our 
own sample, we only choose clusters that are in a contiguous 
region of the survey, and this limits the number of clusters 
in our local comparison sample to 22, at 0.08 < z < 0.09. 

In addition, we repeat the same measurement for the 
rich z = 0.023 cluster, Coma, using the same SDSS data, 
with the red-sequence galaxies selected using (u — g) colour 
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Figure 13. Redshift distribution of clusters in our sample, split 
into two subsamples based on 2Vfedm*+2- The solid histogram 
represents the distribution of the subsample of clusters with 
N re d m* +2 > 20, and the dotted histogram represents those with 
12 < N Tec [ m *+2 < 20. The distribution is smooth with redshift, 
peaking at around z ~ 0.3. 



to bracket the 4000A break at its redshift. Note Coma has 
a iV re d,m*+2 ~ 30 and A^oo ~ 100 (see £|4.1[1 . and therefore 
is ~ 4 times richer than the typical clusters in our sample 
and the Yang07 sample. 



6 LUMINOSITY FUNCTION CONSTRUCTION 

In this paper, we focus on the r' — band red-sequence lumi- 
nosity function (LF) in the core regions of clusters over the 
redshift range of 0.17 < 2 < 0.36. We choose the r' band 
because it has the deepest photometry, and it is a red band 
which is less sensitive to recent star formation than bluer 
bands. We divide our cluster sample into three redshift bins 
and stack all clusters in each bin to obtain a composite red- 
sequence LF, to reduce the noise due to the uncertainty 
in cluster membership determination and cluster-to-cluster 
variation. The width of the redshift bin is chosen in a way 
that the number of clusters in each bin is roughly the same, 
to give similar statistics. 

fc-corrections or (k + e) corrections are applied when 
necessary. T he corrections are calculat ed using the same old, 
single-burst iBruzual fc Chariot] (|2003T l model used to define 
the CMR, and are magnitude dependent. The fc-correction 
is about 0.2 — 0.5 mag at z < 0.27 (the two lower redshift 
bins), and about 0.5 — 0.7 mag at z ~ 0.36 (the highest 
redshift bin). The (k ± e) correction is less than ~ 0.2 mag 
at z < 0.36. 

We construct a composite CMD in each redshift bin to 
define the red-sequence. For clusters in each colour slice we 
calculate the position of every galaxy in the CMD relative 
to m* in that colour slice, and shift those that fall onto the 
red-sequence defined by the wide (r' — i') colour slice (to 



eliminate obvious foreground and background contamina- 
tion) to the central redshift in both colour and magnitude. 
The first column in Figure Q3] shows the composite CMDs 
of galaxies that are within a radius of 0.5 Mpc from cluster 
centres at three different redshifts. The dashed lines indicate 
±0.2 mag from the modelled (g' — r') vs. r' relation. To more 
accurately define the red-sequence, we fit a new (</ — r') vs. 
r' relation using galaxies brighter than m* + 2 (m* is in- 
dicated by the pentagon) that are within ±0.2 mag from 
the model (g' — r') colour (as indicated by the blue points). 
We divide those galaxies into magnitude bins of 0.5 mag, 
calculate the median of the colour distribution in each mag- 
nitude bin, and fit the medians as a function of magnitude 
to a straight line, indicated by the central solid lines. To 
examine the fit more closely, we subtract the fitted (</ — r') 
vs. r' relation and plot the relative position of each galaxy 
to the fitted CMR, which is shown in the second column in 
Figure 1141 with the histograms in the third column show- 
ing the distribution of the residuals (down to m* ± 2) . The 
resulting CMD is centred at zero colour-difference relative 
to the fit by construction. To calculate the width, a, of the 
colour distribution around the fitted (g' — r') vs. r' relation, 
we mirror galaxies that are redder than the fitted CMR to 
avoid contamination from the blue cloud, and apply 3<j clip- 
ping. This a is calculated from galaxies that are brighter 
than m* ± 2, and is not magnitude dependent. The upper 
and lower solid lines in each panel are the 2a bounds of the 
colour distribution of galaxies brighter than m* + 2. 

Isolating red-sequence galaxies is non-trivial; both 
intrinsic scatter and photometric uncertainties can make 
it difficult to cleanly separate the red-sequence and blue 
galaxies. This is particularly a problem at faint magnitudes, 
where photometric errors are large, and the blue cloud 
dominates the population. Here we select red-sequence 
galaxies in several different ways and show how they affect 
the results. The four methods we use are : 

1) Red_4cr, where red-sequence galaxies are defined as those 
that are redder, but not more than 4<r, than the best-fit 
CMR. The total is twice this number. 

2) Red_all, a slight variation of method Red_4cr, where 
red-sequence galaxies are defined as all those that are on 
the red side of the best-fit CMR. Again, the total is twice 
thi s number. Method Red _4q and Red-all ar e mot ivated 
by iBarkhouse et al.l (|2007T ) and Idlbank et all (|2008h (but 
in their work the red-sequence galaxies are not mirrored at 
the bright end). 

3) P10_2<r, where red-sequence galaxies are defined as those 
that have > 10 per cent probability belonging to the colour 
slice defined by the best-fit CMR with a width of ±2er as 
indicated by the solid lines in Figure 1141 This method is 
consistent with how the subs ample in each colour s lice is 
selected (SJ32J prescribed by iGladders fc Yed |2000| . Note 
each galaxy that satisfies this criterion is counted as one 
galaxy, not weighted by this probability. 

4) NP_2er, where red-sequence galaxies are defined as those 
that are completely contained within ±2a from the best-fit 
CMR. 

To correct for the background contamination, the same 
strategy is applied to a sample of galaxies at 500 random po- 
sitions within the four CFHTLS wide fields, for each finely- 
interpolated colour slice. They are shifted to each central 
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redshift bin in the same way as for the cluster+field sam- 
ple, and the average number count is used when subtracting 
the field contribution. The error bars are estimated assum- 
ing that the noise follows Poisson statistics, including the 
background subtraction. Note that the Poisson error on the 
background is negligibly small, because of the large num- 
ber of random fields. We do not include the variance of the 
background field distribution in the error estimate, because 
here we are primarily interested in the average properties of 
a large ensemble of clusters, not the cluster-to-cluster vari- 
ation. 



7 RESULTS 

7.1 Red-sequence Luminosity Functions 

7.1.1 Our sample 

Figure [T5l shows the r'— band red-sequence LFs of all clusters 
with N rf ,d,m'+2 > 12 in the range of projected radius < 
r < 0.5 Mpc, fe-corrected to rest frame, calculated using the 
four methods described in Sjl The depth of the CFHTLS 
data enables us to reach M r / ~ — 17 at z ~ 0.2, a limit that 
has never been probed before at this redshift for such a large 
sample of clusters. 

In all panels in Figure 1151 the dotted histograms are 
the cluster+field counts, the dashed histograms are back- 
ground counts, and the solid histograms with error bars are 
the background-subtracted LFs. The number count in each 
magnitude bin is normalized by the number of clusters con- 
tributing to that bin. The number of clusters in each redshift 
bin is 33, 44 and 50 from low to high redshift, and all clus- 
ters contribute fully to all but the faintest magnitude bin 
due to the completeness limit. 

At bright magnitudes, all four methods give consistent 
results, so it does not matter whether one uses only galaxies 
on the redder side of the CMR or not. However it makes 
a difference at faint magnitudes. In Figure 1151 comparing 
the net counts from the four methods at magnitudes fainter 
than A/ r fc , = -20.5 (2 magnitudes fainter than AT), P10_2cr 
gives the highest count. This is an overestimate because, at 
faint magnitudes, the error on the colour is larger; galax- 
ies that are just outside the colour slice can still have a 
probability of > 10 per cent belonging to the colour slice. 
Note that the LFs are background-subtracted; therefore this 
higher net count in method P10_2<r at faint magnitudes 
is not due to field contamination, but due to the contri- 
bution from the cluster blue cloud. Method NP_2u suffers 
from the same problem as P10_2<7 over the magnitude range 
—20.5 < Mfi < —18.5, as can be seen from the net count 
shown in the bottom row. However, at the very faintest mag- 
nitudes, NP_2<r gives the lowest net count among the four 
methods. This is because, in this case, the width of the fixed 
colour slice is not wide enough compared to the error on 
the colour at this faint magnitude; thus it underestimates 
the number of red-sequence galaxies, due to the preferen- 
tial scattering of the red-sequence galaxies out of the colour 
slice. This problem cannot be alleviated by using a wider 
colour slice because then it would significantly overestimate 
the number of red-sequence galaxies at slightly brighter mag- 
nitudes given the above reasoning about method P10_2<j. 

On the other hand, method Red_4cr and its slight vari- 



ation, Red_all, mirror only galaxies on the redder side of the 
CMR, and thus do not include contributions from the cluster 
blue cloud that occupy regions blueward of the CMD (espe- 
cially at faint magnitudes). Moreover, since method Red_all 
takes account of all galaxies redder than the best-fit CMR, it 
does not underestimate the number of red-sequence galaxies 
regardless of the photometric errors. In method Red_4cr, a 
cut of 4cr redder than the best-fit CMR is applied to elim- 
inate high redshift galaxies with colours much redder than 
the best-fit CMR. This cut corresponds to ~ 0.17-0.27 mag 
from our lowest redshift bin to the highest, which is com- 
parable to or greater than the ler error on the colour even 
at the very faintest magnitude (~ 0.2 mag at r' ~ 23, see 
Figure [2}; therefore, this cut is broad enough to not signif- 
icantly underestimate the number of red-sequence galaxies 
scattered off the CMR due to photometric errors. Compar- 
ing the results from these two methods (first and second 
row in Figure [TSJ), we see that the net counts are consistent 
(solid histograms), but the background (dashed histograms) 
in method Red_4er is significantly reduced compared to that 
in method Red_all; therefore, we conclude that Red_4<r is the 
best, and in the rest of the paper, we will use this method 
for all the analysis. 

The most outstanding feature of our red-sequence LF 
at z ~ 0.2 is a significant and broad dip starting at M r i ~ 
—20.5. The number of red-sequence galaxies reaches its max- 
imum at M r i ~ —20.5, and then decreases to 40 per cent of 
the maximum value at M r i ~ —18.5, and comes back up at 
magnitudes fainter than that. At this redshift, M r i ~ —18.5 
corresponds to r' ~ 22, where the error on the colour and 
total magnitude is still small (see Figure [2}; thus this inflex- 
ion is robust. This feature is also present at z ~ 0.27. It is 
hard to discern at z ~ 0.33 due to the magnitude limit of the 
data and the fact that the error on the colour is getting large 
at the very faintest magnitudes. The dip in the LF proba- 
bly suggests that the LF is made up of a mixture of two 
populations of red-sequence galaxies, possibly giant/regular 
elli ptical and dwar f ellipt icals (dEs): as can be seen in figure 
1 of iBinggeli et all (|l988T l. the LF of E and SO galaxies peak 
at bright magnitudes and drop off at the faint end where the 
LF of dEs starts to increase. The disappearance of dEs over 
the magnitude range seen in our LF could mean that those 
dwarf galaxies on this mass scale are either disrupted (by, 
for example, tidal forces) or they merge into more massive 
galaxies. 

Although the existence of dips in the red- 
sequence /earl y- type LFs has been established in man 



sequence/ early- type L,ts nas Deen es tablished m many 
studies (e.g. Popesso et all 2006|; Barkhouse et all 120071 ; 



iMercurio et al.ll2006l : ISecker fc Harris! 1 19961 ). the depth of 
the dip and the faint-end slope of the LFs are still not well 
constrained, varying significantly in different studies. We 
examine this issure more closely in tj8.ll 

In most cases, instead of fitting a single Schechter func- 
tion, double Schechter functions: one for the bright end and 
one for the faint end; or Gaussian (bright end) + Schechter 
(faint end) functions are used to try to accurately represent 
the shape of the LFs with dips. In our two lowest redshift 
bins, we find that double Schechter functions do not provide 
a good fit to the LF (without removing the brightest cluster 
galaxies). Instead, a Gaussian + Schechter function fit de- 
scribes the shapes of our LFs better. For the highest redshift 
bin, due to the limit of the data, it is hard to tell whether 
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Figure 14. Left panel: Composite CMDs of galaxies that are within a radius of 0.5 Mpc from the cluster centre at three different 
redshifts. The dashed lines indicate ±0.2 mag from the modelled (g' — r') vs. r' relation. Central solid lines are the new best-fit CMR, 
based on galaxies that are brighter than m* + 2 (m* is indicated by the pentagon), shown as the blue points. The upper and lower solid 
lines are the 2a bounds. Right panel: in the left column is the CMD with the best-fit CMR subtracted, and in the right column the 
histograms show the distribution of galaxies relative to the fitted CMR down to m* + 2. The upper and lower solid lines in each panel 
are the 2a bounds. See text for details. 



the inflexion exists; therefore a single Schechter function pro- 
vides an acceptable fit as well. In the lowest redshift bin, all 
six parameters in the Gaussian + Schechter function are 
free; but in the two high redshift bins, the data are not ad- 
equate to constrain the faint end and therefore we fix the 
M* and a at the best-fit values from the lowest redshift bin. 
Note we do not explore the degeneracy among the param- 



eters as that is not our purpose here; we only seek a set 
of parameters that accurately reflect the shape of our LFs. 
The best-fit parameters of our LFs obtained using method 
Red_4<r are listed in Table [3] and the fits are plotted as solid 
curves over the fc-corrected and background-subtracted LFs 
(histograms) in the top panel in Figure [Tol Within the un- 
certainty, the fit to the LF in the lowest redshift bin can also 
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Figure 15. Composite r'— band red-sequence LFs of clusters with N rec i m*+2 > 12 in the projected radius range of < r < 0.5 Mpc, 
calculated using the four methods described in the text, shown in four rows respectively. The top row, Rcd_4<7, is our preferred method, 
as it minimizes the contamination from the blue cloud, and background galaxies. Magnitudes are fc-corrected to rest frame. The dotted 
histograms are the cluster+field counts, the dashed histograms arc background counts and the solid histograms with error bars are the 
background-subtracted LFs. See text for discussions on each method. 
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Figure 16. Top panel: fc-corrected, background-subtracted red-sequence LFs (histograms) obtained using method Red_4cr in three 
redshift bins with the best-fit Gaussian + Schechter functions overplotted as solid curves (dotted curves show the two components 
separately). For the lowest redshift bin, all 6 parameters in the fit are free; for the two higher redshift bins, the M* and a are fixed 
at the best-fit values from the lowest redshift bin. fc-corrected rest-frame absolute magnitude and apparent magnitude are labelled on 
the bo ttom and top axes, respectively. The two vertical dashed lines indicate the division between dwarf and giant as in Lucia et aU 
(2007). Bottom panel: LFs with (k + e)-correction applied. Solid histograms represent (k + e)-corrected LF at each redshift. The one in 
the lowest redshift bin is plotted for reference on top of that in the two higher redshift bins as dotted histograms. Top and bottom axes 
label the apparent magnitude and (k + e)-corrected rest-frame absolute magnitude, respectively. Again, the vertical dashed lines indicate 
the division between dwarf and giant. 
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Table 3. Best-fit parameters of the LFs of our cluster sample con- 
structed using method Red_4t7. Bracketed values are fixed from 
the lowest redshift bin. 





Gaussian 


Schechter 


z 


< M r i > a cf> gau 


M*, a tj> 


0.20 


-20.75 1.3 9.0 


-17.75 -1.15 22.5 


0.27 


-20.75 1.3 8.5 


(-17.75) (-1.15) 26.5 


0.33 


-20.75 1.3 9.0 


(-17.75) (-1.15) 50.0 



fairly represent the LFs in the two higher redshift bins. Note 
the fits are for fc-corrected LFs, but given the low redshift 
probed here, the evolutionary correction is small; thus this 
implies no strong evolution in this redshift range. 

In the bottom panel of Figure [161 we plot our (k + e)- 
corrected LFs. Solid histograms represent the LF in each 
redshift bin, and the one from the lowest redshift is over- 
plotted as dotted histogram in the two high redshift bins. 
There does not seem to be a strong evolution above that 
of passive evolution. In each panel, the two vertical dashed 
lines indicate the div ision between dwarf and giant as in 
iDe Lucia et~ai1 (|2007h . Galaxies brighter than the magni- 
tude indicated by the line on the left are considered giants; 
and galaxies between the two lines are denned as dwarfs (the 
dwarf-to-giant ratio, DGR, will be discussed later on). 



7.1.2 Low-redshift comparison 

As seen in the above section, there does not seem to be a 
strong evolution of the red-sequence LFs within our sample 
at 0.2 <; z <; 0.4. In this section, we compare the LFs of our 
sample with that of the Yang07 sample and the intensely 
studied local rich cluster, Coma. In Figure [T7l the solid his- 
togram is the red-sequence LF of our sample at z ~ 0.2, with 
the brightest cluster galaxies removed (for easier comparison 
with the literature later). Solid points and triangles are the 
red-sequence LFs of the Yang07 sample (z ~ 0.085, J5j> and 
Coma (z = 0.023), respectively, calculated using our method 
Red_4cr, and passively evolved to z = 0.2. To compare the 
shapes of these LFs, in Figure [T7] they are all normalised to 
give the same total number count down to M r i = —20.5. 
Since these LFs are constructed using the same method, 
and the filter combinations used at these redshifts bracket 
the 4000A break in a similar way, the differences in their 
faint ends are not likely due to method differences, but gen- 
uine, suggesting a steepening of the faint end with decreasing 
redshift (but note that Coma is ~ 4 times richer than the 
typical clusters in our and the Yang07 samples, and we will 
discuss this in H8.3.5P . 



7.2 Red-sequence Dwarf-to-Giant Ratio (DGR) 

The DGR is the ratio of the number of faint galaxies to 
bright ones, and thus is commonly used as a simple indicator 
of the shape of the LF. Since red-sequence galaxies represent 
the population that have had their star formation shut off, 
the evolution of red-sequence DGR traces the history of the 
quenching of the star formation. 



Our sample, method Red_4f? (z~0.3) 

• YangO? sample, method lied 4a (7-0.085) 
A Coma, method Red_4<j (z = 0.023) 




Figure 17. r'— band red-sequence LFs at different redshifts, pas- 
sively evolved to z = 0.2. Magnitudes are fc-corrected to rest- 
frame at z = 0. The solid histogram is the red-sequence LF of our 
sample at z ~ 0.2, with the brightest cluster galaxies excluded. 
Solid points and triangles are the red-sequence LF of the Yang07 
sample (z ~ 0.085, jj5} and Coma (z = 0.023), respectively. All 
LFs are calculated using our method Red_4cr, and normalized to 
give the same cumulative number count down to M k , = —20.5. 
The two vertical dashed line s show the appr o ximat e division be- 
tween dwarf and giant as in IDe Lucia etatl ||2007|) . See text for 
details. 



7.2.1 DGR of Our Cluster Sample 

To make it easier to compare with previous works, here 
we k + e correct the apparent magnitude to rest frame at 
z = 0, and adopt the d efinition of dwarf and giant used 
by IDe Lucia etall i|2007l ). i.e. galaxies brighter than —20 
in rest-frame V in the Vega system are considered as giant 
and those between —20 and —18.2 are defined as dwarf. We 
convert the definition of dwarf and giant in V to our rest- 
frame r' magnitude using the (V — r') colour produced by 
the same old, single-burst model described in i|3.1l The tri- 
angles in Figure[TH]are the red-sequence DGRs of our cluster 
sample with N re d >m *+2 > 12, in the radius range < r < 0.5 
Mpc, calculated using method Red_4cr. The bottom and top 
axis labels show the look-back time and redshift respectively. 
The small error bars on our data points are estimated as- 
suming Poisson statistics, without including the variance of 
the background field distribution, because what we are af- 
ter is the average behaviour of the large sample. To get an 
idea of how large the error bar would be on the DGR of a 
single cluster, we calculate the uncertainty on the DGRs for 
a few individual clusters that have DGRs similar to the en- 
semble average, taking into account the noise on the mean 
background and the variance of the distribution of the 500 
random fields. The typical uncertainties are plotted as dot- 
ted error bars on the triangles in Figure 1181 We also find 
a wide cluster-to-cluster variation, and 8 out of 127 of our 
clusters have DGR > 2. We will explore this further in the 
following discussion section. 
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Figure 18. Red-sequence DGRs as a function of look-back time and redshift. Triangles are the DGRs of our cluster sample with 
N re d m *+2 > 12, in the radius range of < r < 0.5 Mpc, calculated using method Red_4cr. The solid error bars indicate uncertainties 
on the sample average DGRs, while dotted error bars represent typical errors on the DGR of a single cluster whose DGR is similar to 
the ensemble average. Open pentagons are the DGRs of a rich subset of our cluster sample with N ret i m *+2 > 18 (discussed in i|8.3.5l l, 
plotted with slight offset in redshift for clarity. The solid squa re and pentagon are th e red-sequence DGR of the Yang07 sample and 
Coma, ca lculated using met hod Red_4cr. The filled circle is from lBarkhouse et al.l J2007T) ■ Crosses are the DGR of our sample, calculated 
using the I Stott et al.l 1(2007) colour cut. See text for more discussion. 
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7.2.2 Evolution with Redshift 

To explore the redshift evolution of the red-sequence DGR, 
we also calculate the DGR for the local comparison sam- 
ple and Coma (*J5|, shown as the square and solid pentagon 
in Figure 1181 The DGRs we have calculated consistently at 
different redshifts show no strong evolution in the redshift 
range 0.2 ^ z & 0.4, but a rapid increase from z ~ 0.2 to 
z ~ 0. Assuming that the bright end of the LF has not 
evolved significantly, our results imply that the number of 
dwarf galaxies has increased by a factor of ~ 3 over the last 
2.5 Gyr. This rapid evolution since z ~ 0.2 is not incon- 
sistent with the predictions from either single burst simple 
stellar population (SSP) or quenched models bv lSmith et al.l 
(2009) through the measurement of the age of faint red 
galaxies in Coma using line indices (see their figure 15), 
although our results agree more with their prediction for 
the outskirt regions instead of the core regions. Despite the 
rapid increase in the number of dwarf galaxies, they do not 
contribute significantly to the growth of the stellar mass on 
the red-sequence, per cluster. Using the stellar mass-to-light 
ratios from the same simple model we used to construct the 
CMD, an increase in the DGR from ~ 0.8 to ~ 2.5 corre- 
sponds to a ~ 15 per cent increase in the stellar mass on the 
red-sequence, per cluster. 



8 DISCUSSION 

8.1 Comparison with Literature 

Studies on cluster red-sequence LF in the literature have 
used quite different colour cuts to select red-sequence galax- 
ies. In this section, we compare our red-sequence LFs and 
DGRs with that in the literature. Given the redshift evo- 
lution shown above, we split the comparison into two red- 
shifts. The LFs in the literature are transformed into our r' 
filter and passively evolve d to z — 0.2, using the same SSP 
iBruzual fc Charlotl <|2003l ) model described in i]3.1l Magni- 
tudes are converted to the cosmology used in this work, if 

different. 

In the work by Bar khouse et al.l ((2007) , the red- 
sequence galaxies are selected in a way that is very similar 
to our method Red_4cr (see Sj6]); therefore it is directly com- 
parable to our measurement of the Yang07 sample (Figure 
I17[) given the similar redshift. In the left panel of Figure [191 
the solid histogram is the red-sequence LF of a sample of 
Abell cluste rs that are X-ray detec ted in t he EINSTEIN X- 
ray im ages (jJones fc Forman|[l999T) used in lBarkhouse et ail 
(2007), measured within the same physical radius (< 0.5 
Mpc) as in this work (Barkhouse, priv. comm.). The filled 
circles represent the LF we measured for the Yang07 sam- 
ple. The two LFs agree fairly well with each other in gen- 
eral. Note the LFs are normalized to give the same cumula- 
tive number count down to M r i = —20.5 (around m* + 2); 
if a slightly fainter magnitude cut is used, given the rel- 
atively large uncertainties at the bright end, the two LFs 
would then have an even better agreement at the faint mag- 
nitudes. A simple way to circumvent the issue of the rela- 
tive normalization is to calculate the DG R. The DGR cal- 
culate d from the LF of the sample from iBarkhouse et all 
i|2007h is shown as the filled circle in Figure 1181 and it is 
consistent with that from the Yang07 sample (the square). 



At a similar redshift, IStott et all (|2007n studied 10 X-ray- 
luminous (Lx > 5 x 10 44 ergs s _1 ) clusters at z = 0.13 
from the Las Campanas/AA T Rich Cluster Survey (LARCS; 
IPimbblet et~ail l200ll . 120061 ). and their LF is shown as the 
open circle s in Figure | 19l It has a s teeper faint end than that 
from both Bar khouse et al.l (|2007l ) and the Yang07 sample. 
We attribute this difference to the way red-seq uence galaxies 
are selected: in the work bv lStott et al.l (|2007T l. red-sequence 
galaxies are considered as those within a fixed colour cut 
around the CMR (\A(B - R)\ < 0.4 mag). This cut is very 
wide, and therefore would include contamination from the 
blue cluster members, overestimating the number of faint 
red-sequence galaxie s. Other work s at a simil ar redshift in- 
cludes , for example, iTanaka et al.l l|2005l ) and lHansen et al.l 
(|2007h . We do no t compare our LFs wi th that in those works 
here because, in ITanaka et al.l l|2005l ) the cluster environ- 
ment is defined using local density, and thus is not clear 
within what p hysical radius the m easurement is made. And 
in the work bv lHansen et al.l (|2007l ) the cluster sample spans 
over a redshift range (0.1 < z < 0.3) within which we see 
evolution in this work; therefore it is not directly compara- 
ble. 

In the right panel of Figure 1191 the solid histogram is 
the LF of our sample at z ~ 0.2, calculated using method 
Red_4<r, with the brightest cluster gal axies included for easy 
comparison with the study from lSmail et all (|l998l ). Circles 
and triangles are red-sequence LFs of a sample of 10 X-ray 
clust ers from ROSAT A ll Sky Survey at z = 0.22 - 0.28 
from Ismail et all fl998). Circles represent red galaxies se- 
lected using (U — B) colour with a width from 0.28 — 0.43 
mag as a function of magnitude, and triangles are selected 
using (B — I) colour with a width from 0.18 — 0.33 mag. 
The filter combination of (B — I) is similar to (g' — r'), and 
the relatively narrow colour slice with varying width as a 
function of magnitude makes it more similar to our method 
P10_2cr. The LF selected using (U — B) has a better agree- 
ment with our LF from method Red_4er, which could be due 
to the fact that at z ~ 0.2 the (U — B) colour is more sen- 
sitive to recent star formation than (B — I), and therefore 
the population that is red in (U — B) is less contaminated 
by the cluster blue cloud. As a result, the DGR from this 
LF would be located close to that of our sample at z ~ 0.2 
if plotted in Figure 1181 To show how a wide colour cut af- 
fects the LF and DGR of our samp le, we calculate th e LF 
of our sample using the colour cut of lStott et al.l (|2007h . We 
transform the difference in (B — R) into (</ — r') making 
use of the colour difference between E and Sab galaxies in 
(B — R) and in (</ — r') at z = 0.2 from iFukugita et al.l 
\ 19951 ). The resulting LF for the sample at z ~ 0.2 is shown 
as the dotted histogram in the right panel of Figure 1191 and 
the DGRs are shown as crosses in Figure 1181 Comparing 
with our LF calculated using method Red_4<r, it produces 
significantly more faint red-sequence galaxies, and as a re- 
sult gives a higher DGR. This shows that the faint end is 
sensitive to the selection method; therefore when looking for 
evolution with the redshift, it is important to make sure all 
methods are consistent. 

The DGR of Coma we calculated (2.5 ± 0.8) is compa- 
rable to that in other studies in the lit erature that used sta- 
tistical background subtraction (e.g. iDe Lucia et al ] 120071 : 
IStott et al]l2007l : lAndreonll2008l ). Also, a complete spectro- 
scopic study on Coma by Marzke (priv. comm.) obtained a 
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red-sequence DGR of ~ 2.1, which confirms the results from 
statistical studies. 

The red-sequence DGRs fro m the studies in the litera- 
ture are compiled in the work bv lGilbank fc Baloghl £2008) 
(see their figure 1). In the redshift range 0.1 < z < 0.4, our 
measurements of the DGRs are systematically lower than 
those in the literature (but note most of those studies do 
not directly overlap with this redshift range). Given that 
the available studies in the literature for the compilation 
all used different selection criteria for red-sequence galaxies, 
and we have shown that the faint end and thus the DGR is 
sensitive to the selection method, we suggest the discrepancy 
is mostly due to how red-sequence galaxies are selected. 

8.2 The Relation between Red Galaxies and 
Passive Galaxies 

Our purpose of studying red-sequence galaxies is to under- 
stand how their star formation is shut off as they move from 
the blue cloud; therefore we want to select a population that 
is as close to a truely quiescent population as possible. 

The population selected using colour that brackets the 
4000A break probably is very close to a population that is 
truely dead, however i t can be contami nated by dusty star- 
forming galaxies (e.g. IWolf et all 120081 ). Since FIR (24/xm) 
emission is sensitive to dusty star formation, it can be used 
to estimate the contamination from the dusty star-forming 
galaxies with red colours, and to select passive populations. 
The CFHTLS Wl field overlaps with the Spitzer Wide-area 
Infrared Extragalactic Survey (SWIRE) XMM-LSS field, 
but the depth of the 24 fim data is not enough to estimate 
the contamination of the dusty star-forming galaxies down 
to low enough SFR at the redshifts prob ed in this work. Ac- 
cording to the work bv lWolf et all l|2008D . in the core of clus- 
ters, the contamination of dusty star-forming galaxies with 
SFR down to 0.14 M, 



yr 



is ~ 20 per cent, dominating at 

10 i 



M* > 10 10 M Q , and < 10 per cent at M* < 10 1U M Q . How- 
ever, the contamination rate is likely to be lower because in 
that work not enough information is available to eliminate 
the contribution from AGNs. Therefore, we believe our LF is 
an accurate representation of the truly inactive population. 

Other than colour, galaxies can be categorized accord- 
ing to their morphology as well. Traditionally, elliptical and 
lenticular galaxies are termed early-type galaxies, while spi- 
ral galaxies are considered late-type galaxies. Also, emission- 
line strength can be used to broadly divide galaxies into 
emission-line (star-forming) or non-emission-line (non-star- 
forming) galaxies. These methods may not select exactly the 
same population. The morphological selection does not al- 
ways trace star formation: although elliptical galaxies are 
generally considered as non-star-forming, observations have 
indicated that a fra ction of elliptical ga laxies do show recent 
star formation fe.g. lKavirai et al]|2007l ); thus they would be 
included in a morphologically selected sample. On the other 
hand, there could be spiral galaxies that do not have star 
formation, as indicated by the partially independent trans- 
formation in morph ology and star forma tion rate in cluster 
environments (e.g. iKodama et al.l [2004t ): thus these spiral 
galaxies would be missed. Therefore, we compare the pas- 
sive pojjulatiori_sj3l^^ these different methods. 

iGodwin fc Peach! (|l977h constructed the LFs of Coma 
for different morphological types using the magnitude mea- 



sured from photographic plates. Based on the numbers 
read off their figure 4, the DGR of the elliptical+lenticular 
galaxies is ~ 2.2. This is consistent with the DGR of the 
red-sequence galaxies in Coma select ed using colour (this 
work, 2.5 ±0.8; Marzke (priv. co mm. ); |Pe Lucia et alj|2007l : 
IStott et ailbOOTl; lAndreonll2008l) 



Christlein fc Zabludofj i[2003l ) constructed R— band LF 



for six Abell clusters at 0.013 < z < 0.067, and split them 
into emission-line and non-emission-line samples, based on 
their [Oil] strength. The DGR we calculate from their LF 
of the non-emission-line cluster galaxies (2.0 ± 0.2), is be- 
tween that of the Yang07 sample and Coma, consistent 
with the trend seen in Figure [18] (although note the crite- 
rion of EW[OII]<5A used to select non-cmission-line galax- 
ies in their work still allows a moderate level of star for- 
mation). These comparsions suggest that at low redshift, 
colour, morphology or emission-line selects mostly the same 
passive p opulation. This is consisten t with the studies by for 
exa mple. ISchweizer fc Seitzerl (1 19921 ); IStrateva et al.l (|200ll ) 
and lHogg et al l |2002i )~ 



8.3 Possible Systematics 

8.3.1 Local vs. Global Background 

We have used a global background subtraction to construct 
the LFs described in £17.11 However, since there are associ- 
ated large-scale structures around clusters, the background 
around clusters might be different from the random field. 
Therefore, we test this effect by constructing the LFs for 
our sample using global background subtraction. Instead of 
taking 500 random fields, for each cluster we take the galaxy 
count in an annulus of radius 10 Mpc around the cluster 
centre, with an area equivalent to the r < 0.5 Mpc core, 
as the background. The resulting LFs are consistent with 
those shown in Figure [151 Thus, we conclude that the use of 
a global background does not bias our results. 

Even if our background is underestimated due to the 
presence of large-scale structure, this will not have a signif- 
icant effect on the red-sequence DGR. As can be seen from 
Figure [T3] using method Red_4cr the background is negligi- 
ble except at low redshift and at magnitudes fainter than 
those used to define the dwarf population. 

Note that the LFs we derived, and those in the literature 
with which we compared, are of the projected galaxy distri- 
bution. Thus, they include galaxies from r > 0.5 Mpc. This 
may lead to a faint end slope which is somewhat steeper than 
that of the genuine core population, but would not signif- 
icant ly affect the DGR measurement (e.g. iBarkhouse et al] 
l2007h . 



8.3.2 Systematic Uncertainties in Photometry 

In this section, we examine any systematic uncertainties in 
photometry. We make use of the data in Dl to examine 
whether the shallower faint end of our LFs is due to incom- 
pleteness, or biased photometry. Two clusters in our sample 
are in both Wl and Dl. We construct the LF based on the 
data from Dl and compare it with that constructed using 
data from Wl. The two LFs are consistent within the un- 
certainty, although the error bars are rather large, given the 
small number of clusters. Thus, the lack of low-luminosity 
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Figure 19. Comparison of our red-sequence LFs with those in the literature. Left p anel: the solid histogram , filled circles, and open 
circles represent the red-sequence LFs of the Yang07 sample. iBarkhouse et all ll2007fl and IStott et al.l ll2007f) , all passively evolved to 
2 = 0.2. Right panel: the solid histogram is the LF of our samp le at z ~ . 2, calc ulated using method Red_4<r, with the brightest 
cluster galaxies in cluded for easy comparison with the study from lSmail et aL I l ll998t) . Circles and triangles are red-sequence LFs from 
ISmail et all dl998h . with red-sequence galaxies selecte d using (U — B), a nd (B — I) colour respectively. The dotted histogram is the LF 
of our sample derived using the colour cut adopted bv lStott et al.l 1120071) . Again, Magni tudes are fc-corrected to rest-frame at z = 0. The 
two vertical dashed lines show the approximate division between dwarf and giant as in |Pe Lucia et al ] (|2007h . See text for discussion. 



galaxies at 0.2 < z < 0.4 we find is not due to insufficient 
depth or biased photometry. This is not unexpected, as we 
have shown we are 100 per cent complete, with small errors 
at the relevant magnitude, based on comparison with the 
same Dl field. 

As an external check, we extract from the SDSS archive 
a subset of galaxies in a 2x2 square degree region that over- 
laps with the CFHTLS Wl field, and compare the total mag- 
nitude (mag_auto for CFHTLS, and mag_model for SDSS) 
of the common galaxies. After transforming the CFHTLS 
magnitude to SDSS magnitude using the formula provided 
on the Terapix websiteTl, the median of the distribution of 
the magnitude difference in g' , r' and i' filters are all less 
than ~ 0.05 mag, and with no obvious trend with magni- 
tude. 

Therefore, the uncertainties in the CFHTLS photome- 
try is not a major source of error. 

8.3.3 Redshift Accuracy 

As described in i]4.3l the redshifts we assign to our clusters 
are accurate to the ~ 0.01 level. If the estimated redshift of 
a cluster is slightly higher or lower than what it actually is, 
the absolute magnitudes of galaxies in that cluster converted 
from apparent magnitudes will be brighter or fainter than 
what they should be, and the intrinsic shape of the LF will 
be distorted. Here we examine how big this effect is using 

t http://terapix.iap. fr/rubrique.php?id_rubrique=241 
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the Schechter function fit from IBarkhouse et"afl (|2007f ). We 
simulate a sample of 100 clusters at z = 0.2, each of them 
with the same intrinsic LF as the Schechter function fit from 
IBarkhouse et al.l (|2007f ). We use a Monte Carlo method, ran- 
domly assigning a redshift to each cluster in the range of 
z = 0.2 ± 0.01. This is equivalent to randomly shifting the 
LF horizontally toward brighter or fainter magnitudes by 
0.29 mag. We then take the average of the 100 shifted LFs as 
the stacked LF for this sample. The stacked LF only makes 
~1.3 per cent difference to the bright end (brighter than 
M T i = —21) of the undistorted LF, and ~0.1 per cent to the 
faint end. Therefore, this does not affect our results. Note, 
this uncertainty of Az ~ 0.01 is comparable to the redshift 
interval between two adjacent colour slices; thus the effect 
of projection f >i8.3.6|l is of comparable uncertainty. 



8.3.4 Aperture Effect 

As mentioned in i]2.2l the colour is measured within a fixed 
aperture of diameter 4.7 arcsec. To examine how colour gra- 
dients of galaxies affect our results, we measure the colour 
in a smaller aperture, of diameter 3 arcsec, and compare the 
colours measured in the two different apertures. The dis- 
tribution of the colour difference of the red-sequence galax- 
ies in the two apertures is symmetric in general. The scat- 
tering of the distribution is smaller than or comparable to 
the uncertainty on the colour determined using the overlap- 
ping pointings (see Figure [2}. At the bright end, the median 
of the colour in the larger aperture is ~ 0.03 mag bluer 
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(smaller than the width of the colour slice). For the LFs 
constructed using methods P10_2<r and NP_2cr, the number 
of red-sequence galaxies might be slightly underestimated 
due to this scattering. However, for methods Red_4cr and 
Red_all, the results are not affected as long as the scattering 
is symmetric. 

8.3.5 Richness Dependence 

Since we detect clusters using only galaxies brighter than 
m*+2, which is around the magnitude that separates giant 
from dwarf galaxies, it is possible that our sample is biased 
towards systems with more giants. To test this, we examine 
whether clusters with fewer giants have higher DGRs. We 
calculate the DGRs for a subset of very rich clusters with 
N re d,m*+2 > 18 (plotted as open pentagons in Figure fTHl 
with slight offset in redshift for clarity), and for a set of 
clusters that are not rich enough to be included in our final 
catalogue (7 < N re d. m *+2 < 10), using our best method 
Red_4cr. We do not find strong dependence of the DGR on 
the richness of the clusters. Therefore, our results are not 
biased due to the detection method. Note the richness of 
our clusters is measured within a radius of 0.5 Mpc from 
cluster centres, and all analysis are carried out for the same 
region. To check if this biases our results, we calculate the 
DGR using a slightly larger radius (r <1 Mpc), and find 
that the change of DGR is within ~ ler of the values shown 
by the triangles in Figure [18] thus it does not change our 
conclusions either. 

For clusters with N re d >m *+2 > 12, we do not de- 
tect obvious correlation between the individual DGR and 
A r re d,m*+2. For the richest cluster in our sample, Abell0362 
at z ~ 0.2 (EH, its red-sequence DGR is ~ 0.7 ± 0.2, 
similar to the sample average DGR. The deficit of dwarf 
galaxies (20.5 < r' < 22.3) is clearly visible on its CMD in 
Figure [5] For the four clusters in our sample that have com- 
parable richness as Coma, their average DGR is 0.9 ±0.1, 
similar to the sample average as well. We also constructed 
the composite red-sequence LF for the four X-ray confirmed 
clusters that have N re d,m'+2 > 12 (but they are all poor 
clusters with kT~ 2 keV), and it is consistent with the LF 
of the whole sample, within the uncertainties. 

Note although Coma is ~ 4 times richer than the typical 
clusters in our sample and the Yang07 sample, the lack of 
strong richness dependence of the red-sequence DGR seen in 
our sample suggests that the rapid evolution between z ~ 
0.2 and z ~ we see in Figure [18] is probably a genuine 
redshift evolution. However, this is based on the assumption 
that the lack of richness dependence still holds at z ~ 0, and 
that Coma is a typical cluster at this redshift, that can be 
fairly compared with our ensemble averages. We will use a 
larger low redshift comparison sample to explore this issue 
in a future paper. 

Note, at higher redshifts, z ~ 0.5, the study by 
iGilbank et all (2008), using a statistical sample of red- 
sequence selected clusters (similar to that presented in this 
work), found that the red-sequence DGR is higher for richer 
systems than t hat for poorer system s. This is contrary to 
the findings of |Pe Lucia et all (|2007T ). who used ~ 10 op- 
tically select ed clusters at z ~ 0.4 — 0.8. The richness 
cuts used in IGilbank et ail (|2008T ) are B gc > 800 for the 
richer subsample (the richness parameter, B gc , is described 



in lYee fc L6pez-CruzHl999T), and 3 00 < B gc < 500 for the 
poorer subsample. De Lucia et al.l (|2007h calculated the dy- 
namical masses (rather than richnesses) via velocity dis- 
persions for individual clusters. The criterion they used to 
separate "richer" (higher velocity dispersion) and "poorer" 
(lower velocity dispersion) clusters was the velocity disper- 
sion being higher or lo wer than 600 km s -1 (which is ap- 
proximately B gc ~ 600. IGilbank et alJliooj ). The division 
used in both studies roughly corresponds to N re d >m *+2 ~ 12 
in our cluster sample, since the virial mass of a cluster with 
a = 600 km s" 1 and r 20 o ~ 1 Mpc is ~ 1 x 10 14 M©. Thus 
the richness/mass-dependence of the red-sequence DGR at 
z <; 0.4 is still an open question. 



8.3.6 Projection Effect 

If two clusters are close enough in redshift space (within 
neighbouring colour slices) , and are projected along the line- 
of-sight, their red-sequences would not be distinguished. As 
a consequence of this projection effect, the number count of 
galaxies per cluster would be overestimated, and the shape 
of the LFs would be slightly smoothed (although this latter 
effect is small, see i|8.3.3[l . Another consequence is that it 
might reduce any trend of the DGR with richness, since a 
rich system that consists of two projected poorer systems 
would give a DGR that is similar to the poorer sys te ms. 

According to, for example, IGilbank et al. I (|2008h . 
for clusters in th e Red-Sequence Cluster Survey 
l|Gladders fc Yedl2005h over the redshift range 0.4 < z < 1, 
the contamination rate from projected systems is ~ 5 — 10 
per cent. Since the method we use to detect our clusters is 
based on the CRS method, we expect the contamination 
rate to be similar. In a later paper, we will use the full 
CFHTLS data to explore this more thoroughly. 



8.3.7 Slope of the CMR 

Since our method Red_4cr mirrors galaxies on the redder 
side of the best-fit CMR to avoid contamination from the 
cluster blue cloud, it is important that the best-fit CMR goes 
through the centre of the red-sequence. Because of the heavy 
contamination from the blue cloud at faint magnitudes, the 
fit to the CMR is carried out at brighter magnitudes down 
to m* + 2, and extrapolated to fainter magnitudes. To test 
the effect of the fitted slope on the DGR, we fix the colour 
of the best-fit CMR at m* , vary the slope by ± ~ 20 per 
cent, and recalculate the red-sequence DGR using method 
Red_4cr. The variation of the DGR due to this slope change is 
about 20 per cent, and thus does not change our conclusions. 

The extrapolation of CMR from bright magnitudes to 
fainter ones implicitly assumes that the slope is independent 
of magnitude. This is reasonable for a truely passive popu- 
lation that is already dead, which is what we are after here. 
For galaxies that are in the transition of moving from the 
blue cloud onto the faint end of the red-sequence, we would 
expect the colour of the red-sequence galaxies to be bluer 
towards fainter magnitudes; therefore a constant slope over 
the whole magnitude range would underestimate this popu- 
lation in transition. Which method to use depends on which 
population one is after. 
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8.3.8 Filter Transformation 

Since the division between dwarf and giant galaxy is de- 
fined in the z = re st-frame V band in the Vega system 
(|De Lucia et al.ll2007t ). we now examine the effect of the un- 
certainty in passband transformation on the red-sequence 
DGR. 

As mentioned in i]3.1l we use a simple single-burst 
model, calibrated with the data of Coma, to produce the 
red-sequence CMD. The (V — r') colour we use to transform 
the definition of dwarf and giant in V to r' is produced by 
the same model. The uncertainty associated with this trans- 
formation could come from a few sources: for example, the 
magnitude-metallicity relation fitted to reproduce the CMR 
of Coma, which, if slightly different, would produce a differ- 
ent colour for a given magnit ude; the uncertain ties on the 
photometry of the Coma data (|Bower et al.| [l992') could also 
introduce an error on the (V — r') colour. However, given 
that the spectral energy distribution of the old galaxy pop- 
ulation is well known, and over the redshift range studied in 
this work the observed r' band is similar to the rest-frame 
V band, we expect the uncertainty on (V — r') to be small, 
at a level of < 0.1 mag. 

We now examine how an uncertainty of ±0.1 mag on 
the division between dwarf and giant in r' would affect 
our results. To maximize the effect, we make the dividing 
point between dwarf and giant 0.1 mag brighter, and the 
faint limit of dwarf 0.1 mag fainter. The resulting DGRs for 
our cluster sample calculated exactly the same way as in 
SOU] are 0.8 ± 0.05 at z ~ 0.2; 1.0 ± 0.05 at z ~ 0.27, and 
1.1 ±0.06 at z ~ 0.33. Note this demonstrates the maximum 
effect by taking an extreme case (two offsets of 0.1 mag in 
directions chosen to maximise changes in DGR); the real 
impact on the results due to the filter transformation is 
smaller than this. 

None of the systematics discussed in this section sig- 
nificantly affects our results. Thus, the lack of evolution at 
0.2 g z % 0.4, and the stronger evolution since z ~ 0.2 seen 
is robust. 



9 CONCLUSIONS 

In this work we detected and compiled a large sample of 
127 clusters over the redshift range of 0.17 < z < 0.36 
from the CFH TLS data using a met hod similar to the 
CRS method (|Gladders fc Yeel [2000), and constructed 
the r 1 — band red-sequence LFs within a projected radius 
of r < 0.5 Mpc from cluster centres. The depth of the 
CFHTLS data enables us to go deeper at these redshifts for 
a larger sample than previous studies. Our main results are 
as follows: 

• We have presented a thorough study of the effect of colour 
selections on the red-sequence LF, and showed that the 
faint end of the LF is very sensitive to how red-sequence 
galaxies are selected. We find that one optimal way to 
minimise the blue cloud contamination is to mirror galaxies 
redder (but no redder than 4er) than the CMR, with a 
second colour cut to further reduce the background. 

• The red-sequence LFs of our sample have a significant 
inflexion centred at M r i ~ —18.5 (especially at z ~ 0.2), 



and thus cannot be described by a single Schechter function, 
suggesting a mixture of two populations. 

• By comparing the red-sequence LFs of our sample with 
that of low redshift samples constructed from SDSS, calcu- 
lated consistently using our best method, we showed there 
is a steepening of the faint end with decreasing redshift 
since z ~ 0.2, but no strong evolution between z ~ 0.2 and 
z ~ 0.4. 

• As a result of the evolution of the LFs, the red-sequence 
DGRs we measured show no significant changes over 
0.2 g z % 0.4, but an increase of a factor of ~ 3 from 
z ~ 0.2 to z ~ (over the last 2.5 Gyr). Also, we do not 
see a strong dependence of the DGR on the richness of the 
clusters within our own sample. 

• We thoroughly checked possible systematics, and showed 
none of them significantly affects our results, and thus our 
results are robust. 

The lack of evolution of the red-sequence DGR over 
0.2 g z g 0.4, and the increase since z ~ 0.2, suggest a rapid 
build-up of the cluster red-sequence in terms of the number 
of dwarf galaxies over the last 2.5 Gyr (the corresponding 
stellar mass on the red-sequence has only inc reased by 15 
per ce nt). Taking the time delay assumed in ISmith et al.l 
(2009) for galaxies to move onto the red-sequence (1 Gyr for 
the single burst SSP models, and 0.5 Gyr for the quenched 
models), our results imply a significant fraction of faint red 
galaxies have been quenched within the last ~ 3 Gyr. 

At higher redshifts than probed in this work, the evo- 
lution of the faint end of the red-s e quenc e is still contro- 
versial, see for exami a le IStott et al.l (|2007T ). 157 Lucia et al.l 
l|2007n . and lAndreonl (|2008l ). Given the different colour se- 
lections employed in those studies, and the sensitivity of the 
faint end to the selection method we have demonstrated, in 
a subsequent paper, we will extend our current study and 
measure the red-sequence LFs consistently out to higher red- 
shifts. With the release of new data from CFHTLS, we will 
have a larger cluster sample, and thus better statistics. 
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